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Abstract 

Weyl semimetals may open a new era in condensed matter physics because they provide the 
first example of Weyl fermions, realize a new topological classification even though the system is 
gapless, exhibit Fermi arc surface states and demonstrate the chiral anomaly and other exotic quan¬ 
tum phenomena. So far, the only known Weyl semimetals are the TaAs class of materials. Here, 
we propose the existence of a tunable Weyl metallic state in Moa:Wi_a;Te 2 via our first-principles 
calculations. We demonstrate that a 2% Mo doping is sufficient to stabilize the Weyl metal state 
not only at low temperatures but also at room temperatures. We show that, within a moderate 
doping regime, the momentum space distance between the Weyl nodes and hence the length of the 
Fermi arcs can be continuously tuned from zero to ~ 3% of the Brillouin zone size via changing Mo 
concentration, thus increasing the topological strength of the system. Our results provide an ex¬ 
perimentally feasible route to realizing Weyl physics in the layered compound Moa;Wi_a;Te 2 , where 
non-saturating magneto-resistance and pressure driven superconductivity have been observed. 
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In 1929, H. Weyl noted that the Dirac eqnation takes a simple form if the mass term is set 
to zero [l|: idi'^ = cp- a'h , with a" being the conventional Pauli matrices. Such a particle, 
the Weyl fermion, is massless but is associated with a dehnite chirality. Weyl fermions may 
be thought of as the basic building blocks for a Dirac fermion. They have played a vital role 
in quantum held theory but they have not been found as fundamental particles in vacuum. 
A Weyl semimetal is a solid state crystal that host Weyl fermions as its low energy quasipar¬ 


ticle excitations 


Q 


16| . Weyl semimetals have attracted intense research interest not only 


because they provide the only known example of a Weyl fermion in nature, but also because 
they can be characterized by a set of topological invariants even though the system is not an 
(topological) insulator. In a Weyl semimetal, a Weyl fermion is associated with an acciden¬ 
tal degeneracy of the band structure. Away from the degeneracy point, the bands disperse 
linearly and the spin texture is chiral, giving rise to a quasiparticle with a two-component 
wavefunction, a hxed chirality and a massless, linear dispersion. Weyl fermions have distinct 
chiralities, either left-handed or right-handed. In a Weyl semimetal crystal, the chiralities of 
the Weyl nodes gives rise to topological charges, which can be understood as monopoles and 
anti-monopoles of Berry flux in momentum space. Remarkably, the topological charges in a 
Weyl semimetal are protected only by the translational invariance of the crystal. The band 
structure degeneracies in Weyl semimetals are uniquely robust against disorder, in contrast 
to the Dirac nodes in graphene, topological insulator and Dirac semimetals which depend 
on additional symmetries beyond the translational symmetry 13Ml7l|. As a result, the Weyl 
fermion carriers are expected to transmit electrical currents effectively. Moreover, the trans¬ 
port properties of Weyl semimetals are predicted to show many exotic phenomena including 
the negative magnetoresistance due to the chiral anomaly known from quantum held theory, 
non-local transport and quantum oscillations where electrons move in real space between 
opposite sides of a sample surface |1CH12|. These novel properties suggest Weyl semimetals 
as a hourishing held of fundamental physics and future technology. The separation of the 
opposite topological charges in momentum space leads to surface state Fermi arcs which 
form an anomalous band structure consisting of open curves that connect the projections of 
opposite topological charges on the boundary of a bulk sample. Without breaking symme¬ 
tries, the only way to destroy the topological Weyl phase is to annihilate Weyl nodes with 
opposite charges by bringing them together in k space. Thus the length of the Fermi arc 
provides a measure of the “topological strength” of a Weyl state. 
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For many years, research on Weyl semimetals has been held back due to the lack of 


experimentally feasible candidate materials. Recently, it was proposed t 


lat 

a 

18 

-21 


a family of 


2l| . Shortly 


after the theoretical prediction, the first Weyl semimetal was experimentally discovered in 
TaAs 1^. So far, the TaAs class of four iso-electronic compounds remains to be only 


experimentally realized Weyl semimetals 


Q,£, 


22 


23|. 


Tungsten ditelluride, WTe 2 , has an inversion symmetry breaking crystal structure, and 
exhibits a compensated semi-metallic ground state 24j-l27l|. The coexistence of inversion 
symmetry breaking and semimetallic transport behavior resembles the properties of TaAs 
and hence suggests a possible Weyl semimetal state. Here, we propose a tunable Weyl 
metallic state in Mo-doped WTe 2 via our hrst-principles calculation, where the length of the 


Fermi arc and hence the topological strength o: 
function of Mo doping. A very recent paper 


the system can be adiabatically tuned as a 
predicted the Weyl state in pure WTe 2 but 


the separation between Weyl nodes was reported to be beyond spectroscopic experimental 
resolution. We demonstrate that a 2% Mo doping is sufficient to stabilize the Weyl metal 
state not only at low temperatures but also at room temperatures. We show that, within a 
moderate doping regime, the momentum space distance between the Weyl nodes and hence 
the length of the Fermi arcs can be continuously tuned from zero to ~ 3% of the Brillouin 
zone size via changing Mo concentration, thus increasing the topological strength of the 
system. Our results present a tunable topological Weyl system, which is not known to be 
possible in the TaAs class of Weyl semimetals. 

We computed the electronic structures using the projector augmented wave method 

le VASP package Ml] within the generalized gradient approx- 


29l. l30| as implemented in t 


imation (GGA) schemes 32j. For WTe 2 , experimental lattice constants were used 33|. 


For MoTe 2 , we assumed that it has the same crystal structure as WTe 2 and calculated the 
lattice constants self-consistently (a = 6.328A, 6 = 3.453A,c = 13.506A). A 8 x 16 x 4 
MonkhorstPack A:-point mesh was used in the computations. The spin-orbit coupling ef¬ 
fects were included in calculations. In order to systematically calculate the surface and 
bulk electronic structure, we constructed a hrst-principles tight-binding model Hamilton for 


both WTe 2 and MoTe 2 , where the ti ght 


orojecting onto the Wannier orbitals 


34 


oinding model matrix elements were calculated by 


36|. which used the VASP2WANNIER90 interface 


37l |. The electronic structure of the samples with hnite dopings was calculated by a linear 
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interpolation of tight-binding model matrix elements of WTe 2 and MoTe 2 . The surface state 
electronic structure was calculated by the surface Green’s function technique, which com¬ 
putes the spectral weight near the surface of a semi-inhnite system. We used W (Mo) s and 
d orbitals and Te p orbitals to construct Wannier functions without using the maximizing 
localization procedure. 

WTe 2 crystalizes in an orthorhombic Bravais lattice, space group Pmn2i (31). In this 
structure, each tungsten layer is sandwiched by two tellurium layers and form strong ionic 
bonds. The left panel of Fig. [11(a) shows a top view of the lattice. It can be seen that the 
tungsten atom is shifted away from the center of the hexagon formed by the tellurium atoms. 
This makes the in-plane lattice constant along the x direction (a) longer than that of along 
the y direction (b). The WTe 2 sandwich stacks along the out of plane z direction, with van 
der Waals bonding between layers (th e right panel of Fig. [11(a)). We used the experimental 
lattice constants reported in Ref. 33], a = 6.282 A, b = 3.496 A, c = 14.07 A. The bulk 
Brillouin zone (BZ) and the (001) surface BZ are shown in Fig. [11(b), where high symmetry 
points are noted. In Fig. [11(c), we show the bulk band structure of WTe 2 along important 
high symmetry directions. Our calculation shows that there is a continuous energy gap near 
the Fermi level, but the conduction and valence bands have a hnite overlap in energy. The 
band gap along the F — R direction is much smaller than that of along the F — W direction 
or r — Z direction, consistent with the fact that the lattice constant b is much smaller than a 
and c. At the Fermi level, our calculation reveals a hole pocket and an electron pocket along 
the F — y direction (in Fig. [11(c)), which agrees with previous calculation and photoemission 
27|. We also calculated the band structure of MoTe 2 by assuming that it is in 


results 


24 


the same crystal structure. As shown in Fig. [11(d), the general trend is that the bands are 
“pushed” closer to the Fermi level. For example, in MoTe 2 , there are bands crossing the 
Fermi level even along the F — A and T — Z directions. We emphasize that, according to 


available literature 


33 


38| , MoTe 2 has a different crystal structure, either hexagonal [2 


or 


monoclinic 33[ both of which have inversion symmetry, but not orthorhombic. Thus in our 
calculation we assumed that MoTe 2 has the orthorhombic crystal structure as WTe 2 and 
obtained the lattice constants and atomic coordinates from hrst-principle calculations. Very 
recently, a paper 39j claimed that MoTe 2 can be grown in the orthorhombic structure. This 
still needs to be further conhrmed. 

We now calculate the band structure of pure WTe 2 throughout the bulk BZ based on the 
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lattice constants reported in [33| . Onr resnlts show that pure WTe 2 has a continuous energy- 
gap throughout the bulk BZ without any Weyl nodes. The k point that corresponds to the 
minimal gap is found to be close to the T — y (Fig. [2](a)) axis. The minimal gap of WTe 2 
is 0.9 meV (Fig. |2](c)). We note that the discrepancy between our results and Ref. j28| is 
due to the sli ghtl y different values of the lattice constants 33|, l40|. The lattice constants 
used in Ref. |28| were at low temperatures Thus the results 28| better refelct the 

groundstate (T = 0) of WTe 2 . We used the lattice constants at room temperatures j33| . 
so our results correspond to tl^ state of WTe 2 at elevated temperatures. The difference 

shows from another angle that WTe 2 is very close to the 


between our results and Ref. 


phase boundary between the Weyl state and the fully gapped state. For many purposes, it 
is favorable to have the Weyl state in a material robust at elevated (room) temperatures. 
Here, we use the room temperature lattice constants for all of our calculations at all Mo 
concentrations. We also note that the very small difference of the lattice constant value 
does not play a role except for undoped or very lightly doped samples x < 2, where the 
separation of the Weyl nodes is beyond experimental resolution anyway. 

We propose Mo doped WTe 2 , Moa;Wi_a;Te 2 , as an experimental feasible platform to realize 
Weyl state in this compound. We have shown that pure WTe 2 is very close to the phase 
transition boundary. Therefore, the k splitting between the Weyl nodes would be beyond 
experimental resolution. On the other hand, another very recent paper proposed a Weyl 
state in pure MoTe 2 MJ, but as shown above the existence of the orthorhombic MoTe 2 needs 
to be conhrmed. By contrast, we show that the moderately Mo doped WTe 2 sample have 
a number of advantages, making it experimentally feasible. First, pure MoTe 2 has many 
irrelevant band crossing the Fermi level along the F — X and T — Z directions, whereas 
the band structure of moderately Mo-doped system is as clean as pure WTe 2 (Figs. [T](c- 
e)). Second, as we will show below, a moderate Mo doping leads to a A: space separation 
of the Weyl nodes that is similarly large as pure MoTe 2 . Therefore, we propose the Mo- 
doped WTe 2 as a better platform for studying Weyl physics. Fig. [2](e) shows the evolution 
of the k space distance between a pair of Weyl nodes as a function of Mo concentration. 
Our calculation shows that a 2% Mo doping is sufficient to stabilize the system in the Weyl 
metal state. Also, the distance between the Weyl nodes increases rapidly at the small doping 
regime. At a moderate doping x = 20%, the k space distance is found to be as large as 
0.03^. As one further increases the doping concentration, the distance seems saturated. 
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The distance is about 0.04^ at x = 40%. The energy difference between the pair of Weyl 
nodes is shown in Fig. Et^g). In Fig. Et^d) we show the dispersion along the momentum space 
cut that goes through the direct pair of Weyl nodes as defined in Fig. E^b). It can be seen 
clearly that two singly generate bands, b2 and b3, cross each other and form the two Weyl 
nodes with opposite chiralities. We name the Weyl node at lower energy as W1 and the 
Weyl node at higher energy as W2. Another useful quantity is the energy difference between 
the extrema of these two bands. This characterizes the magnitude of the band inversion, as 
shown in Fig. El^f). It is interesting to note that, in contrast to the k space distance between 
the Weyl nodes (Fig. E](e)), the energy difference between the Weyl nodes (Fig. E](g)) and 
the band inversion energy (Fig. El(f)) does not show signs of saturation as one increases the 
Mo concentration x up to 40%. In Fig. E](h), we show a schematic for the distribution of the 
Weyl nodes in Mo doped WTe 2 . We observe a pair of Weyl nodes in each quadrant of the 
kz = ^ plane. Thus in total there are 4 pairs of Weyl nodes on the kz = ^ plane. 

A critical signature of a Weyl semimetal/metal is the existence of Fermi arc surface states. 
We present calculations of the (001) surface states in Fig. El We choose the 20% Mo-doped 
system, Moo. 2 Wo. 8 Te 2 . FigureEl(a) shows the surface energy dispersion along the momentum 
space cut that goes through the direct pair of Weyl nodes, Wl(-) and W2(-|-), which arises 
from a single band inversion. Our calculation (Fig. El(a)) clearly shows the topological 
Fermi arc surface state, which connects the direct pair of Weyl nodes. The Fermi arc is 
found to terminates directly onto the projected Weyl nodes. In addition, we also observe a 
normal surface state, which avoids the Weyl node and merges into the bulk band continuum. 
Because the W1 and W2 Weyl nodes have different energies, and because W1 is a type II 
Weyl cone 2^, constant energy maps always have finite Fermi surfaces. Hence, visualizing 
Fermi arc connectivity in constant energy maps is not straightforward. Instead of a constant 
energy map, it is possible to use a varying-energy k^^ky map, i.e. Energy= E{kx,ky), so 
that there are no bulk states on this varying-energy map at all k^, ky points except the Weyl 
nodes. Figure El(f) shows the calculated surface and bulk electronic structure on such a 
varying-energy kx, ky map in the vicinity of a pair of Weyl nodes. A Fermi arc that connects 
the pair of Weyl nodes can be clearly seen. We study the effect of surface perturbations. The 
existence of Weyl nodes and Fermi arcs are guaranteed by the system’s topology whereas 
the details of the surface states can change under surface perturbations. In order to do so, 
we change the surface on-site potentials of the system. Physically, the surface potentials can 



be changed by surface deposition or applying an electric field on the surface. Figure [3l^b) 
shows the surface band structure with the surface on-site energy increased by 0.02 eV. We 
hnd that the normal surface state moves further away from the Weyl nodes whereas the 
topological Fermi arc does not change signihcantly. Figure |3]^c) shows the surface band 
structure with the surface on-site energy decreased by 0.11 eV. The normal surface states 
disappear. The Fermi arc also changes signihcantly. Instead of directly connecting the two 
Weyl nodes in Figs. [3]^a and b), a surface state stems from each Weyl node and disperses 
outside the window. We note that the surface states in Fig. [3]^c) are still topological and are 
still arcs because they terminate directly onto the projected Weyl nodes. We illustrate the 
two types of Fermi arc connectivity in Figs, [^d and e). Figure |3]^d) corresponds to the case 
in Figs. El^a and b), where a Fermi arc directly connects the pair of Weyl nodes in a quadrant. 
Figure [3l|^e) corresponds to Fig. Et^d). In this case, Fermi arcs connect Weyl nodes in two 
different quadrants across the Y — Gamma — Y line. The nontrivial topology in a Weyl 
semimetal requires that there must be Fermi arc(s) terminating onto each projected Weyl 
node with a nonzero projected chiral charge and that the number of Fermi arcs associated 
with a projected Weyl node must equal its projected chiral charge. On the other hand, 
the pattern of connectivity can vary depending on details of the surface. The observed 
different Fermi arc connectivity patterns as a function of surface on-site potential provide 
an explicit example of both the constraints imposed and the freedoms allowed to the Fermi 
arc electronic structure by the nontrivial topology in a Weyl semimetal. 

We further study the surface states via bulk boundary correspondence. We note that 
except Figs. Et^b-e), all other hgures corresponds to the case without additional changes 
to the onsite energy. Specihcally, we choose a closed loop in ky) space as shown in 
Fig. El^j). As we mentioned above, the conduction and valence bands only touch at the eight 
Weyl nodes. Thus as long as the loop chosen does not go through these Weyl nodes, there 
is a continuous bulk energy gap along the loop. In the bulk BZ, the chosen rectangular 
loop corresponds to a rectangular pipe along the kz direction. Then topological band theory 
requires that the net chiral charge of the Weyl nodes that are enclosed by the pipe equals 
the Chern number on this manifold, which further equals the net number of chiral edge 
modes along the loop. For example, the rectangular loop a — /3 — 7 — 5 in Fig. [3]^j) encloses 
a Wl(-) and a W2(-|-), which leads to a net chiral charge zero. The energy dispersion along 
this rectangular loop is shown in Fig. El^g). It can be seen that the bands are fully gapped 


9 


without any surface states along /3 — ^ — 6 — a. Along a — /3, there are two surface states, 
SSI and SS2, both of which connect the band gap. Interestingly, we note that these two 
surface bands are counter-propagating although they seem to have the same sign of Fermi 
velocity. This is because the continuous energy gap a — /3 is highly “tilted”. If we “tilt” 
the energy gap back to being horizontal, then it can be clearly seen that the two surface 
bands are counter-propagating, which means that the net number of chiral edge modes is 
zero. Similarly, we can choose other loops. For example, we choose another rectangular loop 
a' — /3' — 7 — h that encloses only the W2(-|-) Weyl node. Because there are no surface states 
along the two horizontal edges and the vertical edge to the right, we only need to study the 
vertical edge to the left, that is the a' — jS'. The enclosed net chiral charge is -|-1, which 
should equal the net number of chiral edge mode along a' — j3'. The band structure along 
this line is shown in Fig. |3l^h). We see that while the surface band SSI still connects across 
the band gap, SS2 starts from and ends at the conduction bands. Therefore, SSI contributes 
one net chiral edge mode whereas SS2 contributes zero net chiral edge mode. Hence there is 
one net chiral edge mode along the this rectangular loop. By the same token, we can choose 
the loop a" — /3" — j — 6, which does not enclose any Weyl node. Consistently, as shown in 
Fig. |3l^i), along a” — (3”, SSI does not appear along this line and SS2 does not connect across 
the band gap. Hence the net number of chiral edge mode is also zero along this rectangular 
loop. 

we study the constant energy contours of the surface states. We emphasize that (1) there 
is a signihcant energy offset between the W1 and W2 Weyl nodes, and that (2) the W1 
Weyl cones are the type H Weyl cone 28(], which means that at the Weyl node energy, 
its constant energy contour consists of an electron and a hole pockets touching at a point, 
the Weyl node. These two properties are very different from the ideal picture, where all 
Weyl cones are normal rather than type H and their nodes are all at the same energy. We 
show below that these two properties make the surface states’ constant energy contours 
quite different from what one would expect naively. Fig. IHd) shows the calculated constant 
energy contour within the top half of the surface BZ at energy Ei, which is between the 
W1 and W2 nodes in energy. We see three bulk pockets. The corresponding schematic is 
shown in Fig. HJ^a). Specifically, we see a big pocket closer to the Y point, which encloses 
two W1 Weyl nodes with opposite chiral charges. We also see two separate small pockets 
closer to the F, each of which encloses a W2 Weyl node. As for the surface states, from 
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Figs, m^d and e) we see a surface state band that connects the two small pockets, each of 
which encloses a W2 Weyl node. This is quite counter-intuitive because we know that the 
Fermi arc connects the direct pair of Weyl nodes, namely a Wl(-) and a W2(-|-) or vice 
versa. We show that there is no discrepancy. Specihcally, we show that the surface band 
seen in the constant energy contours is exactly the Fermi arc that connects the Wl(-) and 
the W2(-|-) Weyl nodes seen in Fig. HJ^c). To do so, we consider the constant energy contours 
at two different energies, Ei and E 2 . According to the energy dispersion (Fig. 111(c)), we see 
that the big bulk pocket in the constant energy map is electron-like while the two small 
bulk pockets are hole-like. Thus as we increase the energy from Ei to i? 2 , the big pocket 
should expand whereas the two small pockets should shrink, as shown in Fig. 111(b). The 
surface state band keeps connecting the two small pockets as one changes the energy. This 
evolution is shown by real calculations in Figs. lll(e and f). The orange line in Fig. 111(b) 
connects the Wl(-) and W2(-|-) Weyl nodes. At each energy, the surface state band crosses 
the orange line at a specihc k point. By picking up the crossing points at different energies, 
we can reconstruct the Fermi arc that connect the Wl(-) and W2(-|-) Weyl nodes shown 
in Fig. 111(c). Therefore, from our systematic studies above, we show that the Fermi arc 
connectivity means the pattern in which the surface state connects the Weyl nodes. This is 
dehned on a varying-energy (A:^,, ky) map where the chosen E(kx, ky) map crosses the bulk 
bands only at the Weyl nodes. If there is no signihcant energy offset between Weyl nodes 
and if all Weyl cones are normal rather than type II, then the connectivity can also be seen 
in a constant energy contour. However, in our case here, one needs to be careful with the 
simplihed ideal picture, that is to study the Fermi arc connectivity from the constant energy 
contour. Because of the energy offset between the Weyl nodes and because of the existence 
of type II Weyl cones, how surface bands connect different bulk pockets in a constant energy 
contour does not straightforwardly show the Fermi arc connectivity. 

Finally, we discuss the tunability of the length of the Fermi arcs as a function of Mo 
concentration x in our Moa;Wi_a;Te 2 system. (1) The undoped x = 0 sample is fully gapped 
according to our calculations (Fig. 111(g)). (2) A very small Mo concentration (~ 0.5%) will 
drive the system to the critical point, where the conduction and valence bands just touch 
each other Fig. 111(h). The length of the Fermi arc is zero, and hence the system is at the 
critical point. (3) As one further increases the Mo concentration x, the touching point splits 
into a pair of Weyl nodes with opposite chiralities (Fig. lll(i)). The Weyl nodes are connected 
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by a Fermi arc. A way to gap the system without breaking any symmetry is to annihilate 
pairs of Weyl nodes with opposite chiralities. In order to do so, one needs to overcome the 
momentum space separation between the Weyl nodes to bring them together in k space. 
For a sample with a given Mo concentration x, this may be achieved by applying external 
hydrostatic pressure. Thus the length of the Fermi arc provides a measure of the system’s 
topological strength. Such a tunability is not known in the TaAs class of Weyl system. This 
highlights the uniqueness of our work. 
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FIG. 1: Crystal structure and band structure, a, A top view and a side view of WTe 2 
lattice. The silver and yellow balls represent W and Te atoms, respectively. WTe 2 crystalizes in 
an orthorhombic Bravais lattice, space group Pmn2i (31). The lattice constants are a = 6.282 A, 
b = 3.496 A, c = 14.07 A according to a previous X-ray diffraction (XRD) measurement 33|. b, 


The bulk and (001) surface Brillouin zone of WTe 2 . c, Bulk band structure of WTe 2 . d, Bulk 
band structure of MoTe 2 by assuming that it has the same crystal structure as WTe2. Note that, 
in fact, according to available literature MoTe 2 has two possible structures 2, both of which 


are different from the crystal structure of WTe 2 . e, Bulk band structure of Moo, 2 Wo, 8 Te 2 . 
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FIG. 2: Tunable Weyl node separation (topological strength) in Mo-doped WTe 2 . a, 

The k space location that corresponds to the minimal energy gap of the pure (undoped) WTe 2 . 
b, The k space locations of the pair of Weyl nodes within the first quadrant of the = 0 plane 
for Moo, 2 Wo, 8 Te 2 . White and Black dots show the negative and positive chiral charges {Cw), 
respectively, c, Band structures of WTe 2 along momentum space cut as defined in pane (a). 
The four bands appeared in this panel are labeled as bl-b4, respectively, d, Band structures of 
Moo, 2 Wo, 8 Te 2 along momentum space cut as defined in pane (b). We find an inversion between the 
b 2 and bs bands, giving rise to the Weyl nodes W1 and W2 with opposite chiral charges Cw e, 
The k space separation between the pair of Weyl nodes W1 and W2 as a function of Mo doping x. 
f, The energy difference between the extrema of the b2 and b3 bands as a function of Mo doping 
X. This characterizes the magnitude of the band inversion, g, The energy offset between the pair 
of Weyl nodes W1 and W2 as a function of Mo doping x. h, A schematic for the distribution of 
the Weyl nodes within a bulk BZ. All nodes are located on the = 0 plane. 
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FIG. 3: Tunable Fermi arc length and inter-connectivity in Moa;Wi_a;Te 2 . 
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FIG. 3: a, Surface and bulk band structure of Moo, 2 Wo, 8 Te 2 along the momentum space cut that 
goes through a direct panel of Weyl nodes, Wl(-) and W2(+), as defined in Fig. EJ^b). A topological 
Fermi arc surface state connects the Weyl nodes W1 and W2. A normal surface state avoids the 
Weyl nodes and merges into the bulk band continuum, b, Same as panel (a) but with the surface 
on-site energy increased by 0.02 eV. c, Same as (a) but with the surface on-site energy decreased 
by 0.11 eV. d, Schematic of Fermi arc connectivity pattern for panels (a and d). e, Schematic of 
Fermi arc connectivity pattern for panel (c). f, Surface band structure on a varying-energy {kx, ky) 
map. Specifically, we choose a different energy on each [kx, ky) point, so that the there are no bulk 
states at all kx,ky except at the locations of the eight Weyl nodes. This is possible because we 
know that the conduction ad valence bands only cross at the eight Weyl nodes, bf g. Surface and 
bulk band structure along the fc—space trajectory a — P — ■j — S defined in panel (j). h,i, Surface 
and bulk band structure along the fc—space trajectories a' — j3' and a' — 13' defined in panel (j). j, 
Schematic of a quadrant of the surface BZ, showing the /c—space trajectories used in panels (g-i). 
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FIG. 4: Iso-energy contour interconnectivity on the surface of Moa;Wi_ 3 ;Te 2 . a, 

Schematic illustration of the surface and bulk electronic structure on a constant energy Ei. The 
shaded areas represent the projected bulk bands whereas the red line show the surface states. 
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FIG. 4: b, Schematic illustration of the surface and bulk electronic structure on two energies Ei 
and E 2 . The band structure on Ei {E 2 ) are shown by the solid (dotted) lines. The energies Ei and 
E 2 are defined in panel (c). c, Band structure along a k space cut that goes through the direct pair 
of Weyl nodes, Wl(-) and W2(+). The dotted lines denote the energies Ei and E 2 with respect to 
the W1 and W2 Weyl nodes, d, Calculated surface and bulk electronic structure on at the energy 
El over the top half of the surface BZ. e, A Zoomed-in view of panel (d) for the area highlighted 
by the blue box. f, The same as panel (e) but at the energy E 2 . g-i, Schematic band diagram to 
show the evolution of the Moa;Wi_a;Te 2 system as a function of Mo concentration x. j, The length 
of the Fermi arc as a function of Mo concentration x. The arc length equals the k space separation 
of the Weyl nodes. 



